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Data from the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) mission 
show Comet 103P/Hartley 2 is a bi-lobed, elongated, nearly axially symmetric comet 2.33 km in length. 
Surface features are primarily small mounds <40 m across, irregularly-shaped smooth areas on the two 
lobes, and a smooth but variegated region forming a “waist” between the two lobes. Assuming parts of 
the comet body approach the shape of an equipotential surface, the mean density of Hartley 2 is modeled 
to be 200-400 kg m -3 . Such a mean density suggests mass loss per orbit of >1%. The shape may be the 
evolutionary product of insolation, sublimation, and temporary deposition of materials controlled by 
the object’s complex rotation. 

© 2012 Elsevier Inc. All rights reserved. 


1. Introduction 

After carrying out an impact experiment on Comet 9P/Tempel 1 
(A’Hearn et al., 2005a) the Deep Impact mission’s flyby spacecraft 
was retargeted to Comet 103P/Hartley 2 as part of an extended 
mission named EPOXI. Compared to Deep Impact’s first target, 
Hartley 2’s nucleus is one-fifth the size yet produces a similar 
amount of outgassing at perihelion (Groussin et al., 2004; Lisse 
et al., 2009; Meech et al., 2011). This high level of activity made 
103P/Hartley 2 an especially interesting target, and for expecta- 
tions that its surface features would differ from those of Tempel 1. 

Closest approach to Hartley 2 was 694 km at 13:59:47.31 UTC 
on 4 November 2010, 1 week after perihelion passage and 
1.064 AU from the Sun. The images confirmed that its surface is 
very different from that of Tempel 1, and indeed from the other 
comets so far visited by spacecraft (A’Hearn et al., 201 1 ; Schultz 
etal.,2011). 
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We use the image data to produce a shape model (refined from 
that reported in A’Hearn et al. (2011)) and surface feature maps. 
These products are then used to investigate the processes effective 
at the comet’s surface, its mean density, and some aspects of its 
evolution. 

2. Data and methods 

Because the Hartley 2 encounter has been described in depth in 
A’Hearn et al. (201 1 ), we summarize only the most pertinent imag- 
ing instrumental details here. The Deep Impact spacecraft and its 
instruments are described in A’Hearn et al. (2005b), Hampton 
et al. (2005), and Klaasen et al. (2008). Most of the work reported 
here has used the EPOXI Medium Resolution Instrument (MRI) 
camera data; these have a pixel scale of ~7 m at the closest ap- 
proach of 694 km. The EPOXI High Resolution Instrument (HRI) 
nominally produces images at ~5 times better resolution than 
MRI, but was found to be out of focus after launch (Klaasen et al., 
2008). Deconvolution (Lindler et al., 2012) can render some of 
these images suitable for studying forms smaller than are visible 
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Table 1 

Characteristics of Hartley 2. 

Area: 5.24 km 2 
Volume: 0.809 ± 0.077 km 3 
Mean radius: 0.58 ± 0.018 km 
Diameter range: 0.69-2.33 km 

Gravity: 0.0019-0.0044 cm s -2 (if mean density = 300 kg m -3 ) 

Model moments: 5.302 x 10 4 , 31.26 x 10 4 , 31.92 x 10 4 m 2 
Model moment ratios: A/C: 0.166 ± 0.004, B/C: 0.979 ± 0.002 
Body model moment orientations: A: 89.73°, 207.56°E; B: 0.05°, 106.14°E; C: 
0.25°, 16.14°E 

Long axis at closest approach: RA: 226.12°, Dec: 39.37° 2010-2011- 
04T13:59:47.7 

Minimum moment at closest approach: RA: 225.92 ± 0.9°, Dec: 39.60 ± 0.7° 
Intermediate moment at closest approach: RA: 330.0, Dec: 16.3, ±12° along 
arc perpendicular to minimum moment 


in the MRI data, but this processing often introduces artifacts that 
complicate their use in geological interpretation. 

Most useful MRI data were obtained within a range of 1600 km 
of the comet (16 m/pixel) and span phase angles of 80-94°. A few 
HRI images were useful out to ranges of ~2500 km. Navigation 
data in the form of “SPICE” kernels (Semenov et al., 2004; Semenov 
and Acton, 2006) are the basis for all geometric work on the comet. 
An important characteristic of this comet is its complex rotation 
(Belton et al., 2012) which limits any mapping convention that 
would be tied to the surface. This rotation includes a primary per- 
iod, increasing through encounter, and a roll about the long axis, 
with a decreasing period. Most useful mapping data were obtained 
within a period when a simplified model of rotation predicts that 
the object orientation changes by ~0.6°. Because of our desired 
mapping conventions (see below), the routine mapping was done 
with a simple rotation model that is likely to introduce relative po- 
sition errors of 0.3° between extremities of the object. This error is 
less than 1 pixel at the ends of the object. 

Determination of the shape and accurate relative positioning of 
the images relies upon stereo control points with image pointing 
adjustments (Thomas et al., 2002). There are 244 manually mea- 
sured points in the ~50% of the comet that is directly observed. 
Good control of much of the rest of the shape is provided by pro- 


nounced nightside silhouettes arising from light scattered from 
the coma. The mean control point residuals of 0.29 pixel, or 
~3 m suggest that mean dimensions of the comet for the ~50% 
constrained by the stereo data are good to ±10 m. Much of the rest 
of the surface determined from silhouettes is good to ±30 m; most 
of the uncertainty arises from the ambiguity of where a limb pixel 
is along the line of sight. Combining the areas and uncertainties, 
we estimate the mean radius (radius of sphere of equivalent vol- 
ume) is measured to 18 m accuracy (1 sigma); see Table 1. 

Mapping of features and projection of images uses image cubes 
that store latitude, longitude, radius, incidence, and emission at 
each pixel. Line and sample coordinates of features marked in 
the images or the image data themselves can then be arbitrarily 
projected. 

Quantities such as gravity and slope are calculated with an as- 
sumed mean density and rotation vectors. We discuss the determi- 
nation of the nucleus’ mean density in detail in Section 5. 
“Topography” is calculated as dynamic height, the surface poten- 
tial energy divided by an average acceleration (Vanicek and Krak- 
iwsky, 1986; Thomas, 1993) which is close to the height above a 
reference equipotential. The dynamic heights also depend upon 
mean density and rotation rate and are studied in Section 5. 

Mapping uses a prominent feature as a prime meridian: a 
dark spot, actually including part of a shaded side and shadow 
cast from a 20-m high mound. The map coordinate system used 
here is oriented on the basis of the bi-lobed shape: the long axis 
is taken to be the latitude direction. Given the complex and chang- 
ing rotation (Belton et al., 2012) latitudes and longitudes cannot be 
tied to any rotation element. Our maps are not intended as IAU 
proposals. 

The cylindrical mapping of such an elongate, bi-lobate object 
has distorted areas no matter the orientation. Here we use both 
simple cylindrical and a pseudo-equal area projection to better 
handle high latitude data. This method simply measures the dis- 
tance along latitude lines from a reference longitude, and scales 
the north-south coordinate to the average length of a longitude 
line. This scheme is particularly suited to a nearly axially-symmet- 
ric object such as Hartley 2 because the length of the longitude 
lines varies only slightly. 



Fig. 1. Stereo view of Hartley 2. Images are mv03421 49303 J5000002 and mv0342149307_5004051 obtained from ranges of 695 and 697 km. Smaller lobe is the northern 
one; Sun is from slightly left of up. Right side is the evening terminator. Zero longitude goes through the dark spot (partly a shadow) mid way up the body, close to the limb at 
the left side. 
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Fig. 2. Shape data and shape, (a) Limb and control point data locations, simple 
cylindrical projection. Large dots are control point solutions, small dots are 
calculated limb and silhouette positions, (b) Shape model colored by relative 
dynamic heights. (For interpretation of the references to color in this figure legend, 
the reader is referred to the web version of this article.) 

3. Overall results 

Global characteristics of Hartley 2 are listed in Table 1. Preli- 
minary values, close to these, were presented in A’Hearn et al. 
(2011); updates to the shape have largely been local changes of 
<20 m from addition of control points and more limb control. 
Fig. 1 presents a stereo view at the best MRI resolution which cov- 
ers most of the visible object. Fig. 2 shows coverage for data that 
control the calculation of the shape, as well as views of the result- 
ing shape model shaded by the dynamic heights. The mean radius 
(radius of a sphere of equivalent volume) of 0.580 ± 0.018 km 
found here is indistinguishable from that calculated before 
encounter by Lisse et al. (2009) using Spitzer data to determine 
the effective radius of the cross-section at the time of observation. 
This bi-lobed object has near rotational symmetry: calculated mo- 
ments, assuming a homogeneous interior, have relative values of 
1.000, 0.979 ±0.002, 0.166 ±0.004. These model moment ratios 
are unusually but not uniquely symmetric. For comparison, the 
model moment ratio B/C for Tempel 1 is 0.931, it is 0.982 for aster- 
oid Eros, and it is 0.980 for satellite Epimetheus. Several other 
Saturnian small satellites’ B/C range from 0.819 to 0.886 (These 
values calculated by PT). That Hartley 2 is not exactly symmetric 
is noticeable in Fig. 1 when carefully viewing the larger lobe. None- 
theless, this object does not have the asymmetries that are seen on 
many impact-shaped objects (Thomas, 1989). Note that an object 
such as Eros can have moment ratios of nearly 1 while having a dis- 
tinctly banana-like shape. 

4. Terrain types and surface features 

The surface of Hartley 2 is approximately divisible into two 
rough lobes connected by a smooth “waist” (Figs. 1 and 2). Surface 
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Fig. 3. Cylindrical projection maps, (a) Mapped features. Note that there are 
separate bright and dark mound symbols; these are qualitative classifications only. 
Darker and lighter portions of the waist are differentiated. Other smooth areas are 
mapped separately, (b) HRI image mosaic, (c) MRI images. 

maps are shown in Figs. 3 and 4. In detail the surface combines 
varying distributions of mounds generally less than 30 m in any 
dimension, and substantially smoother terrain (at 7 m/pixel). 
While the latitudinally confined waist is the largest smooth area, 
other smooth areas of irregular outline occur in parts of the larger 
lobe. 

4.1. Mounds 

The roughness of Hartley 2 derives from the presence of 
mounds, mostly only a few pixels across, and most of which are 
crudely equidimensional. However, a few of these are unusually 
elongated, either nearly normal to the surface or at low angles to 
the surface (Fig. 1). There are a few spires greater than 40 m in 
height, one of which appears in stereo viewing to be angular in 
shape with a probable overhang (upper part of Fig. 1 near termina- 
tor). Views such as Fig. 1 suggest there are albedo differences 
among the mounds, but estimating reliable albedos from margin- 
ally resolved features with a variety of photometric geometries de- 
mands extreme caution. The qualitative mapping did divide 
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Fig. 4. Pseudo equal-area projections of Hartley 2. Left: Mapped features, same scheme as Fig. 3. Right: Image map, combining both HRI and MRI data. 



Size, m 


Fig. 5. Histogram of mound sizes. Data are mean horizontal dimensions; 5 m bin 
sizes. Fall-off below 20 m is probably largely resolution-related. 

mounds by apparent relative brightness (Fig. 3), though for most 
purposes one should lump them together as one population of fea- 
tures. The size distribution of mounds (Fig. 5) suggests a large pop- 
ulation smaller than the available resolution. 


4.2. The waist 

The distinctive waist (Figs. 1 and 6) is visible in images between 
longitudes 90-280°E and latitudes 2°S-22°N to ~60°N. Silhouettes 
in both HRI and MRI data show this shape occurs over nearly the 
full circumference of the object (Fig. 6). The surface is variegated 
at horizontal scales of 15-30 m. The markings in the waist, mostly 
only 2-3 pixels across, are repeatable image-to-image and are 
well-defined in stereo pairs (Fig. 7a). Thus these subtle forms are 
real and not simply noise. Patterns of different albedo also exist 
on the 100-200 m scale and suggest digitate outlines (Fig. 7a), 
especially near the margins. A darker region between ~10 and 
40°N occupies most of the visible longitude range but fades at 
about 70°E. This limitation may reflect an incidence angle effect 
at values >60°. 



z, km 


Fig. 6. Limb scans around the waist. Abscissa is the distance parallel to projection of 
z axis in the image; ordinate is limb position relative to z axis in image, each is offset 
to show relative topography. Three digit numbers are last digits of the mission 
elapsed time portion of MRI image numbers (for example: mv0342149372). 
Vertical dimension exaggerated by 1 .6. 


The smoothness of the waist may be important in questions of 
its formation mechanism. Fig. 6 shows limb profiles measured to 
~0.15 pixels precision (~1 m in the best images) by software spe- 
cifically designed for edge detection (Dermott and Thomas, 1988). 
Repeatable deviations of ~2 m occur within these profiles. As with 
any limb profile the observed edge is a skyline smoother than a 
straight topographic profile over the position of the limb. Most 
deviations are under 1 m. The very high incidence angle region 
near the east end of the waist is limited in the applicability of shad- 
ows to measuring topography because some of this area is shad- 
owed by protuberances that are part of the small lobe. The solar 
direction does have a westward component near the terminator, 
so some of the waist is seen at incidence angles above 86°. Irregu- 
larities in the terminator of ~10-15m suggest topography of at 
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least ~1 m in height in this region. Thus, while the waist is 
smoother than the rest of the object, formative processes must 
cause, expose, or preserve topography of up to 2 m in height 
(and possibly more) on local scales up to 30 m in horizontal extent. 

4.3. Other smooth areas 

The large lobe has several (7 or more) irregularly-shaped areas 
of smoother and slightly darker appearance than the average sur- 
face of Hartley 2 (Figs. 2 and 7b). These darker areas are flanked 
by slightly brighter materials of varying extent. In combination 
with the mounds these areas comprise a sequence from the center 
outward of dark smooth material, brighter smooth areas, and 
mound-dense regions. The stereo pairs (Fig. 7b) show that these 
dark smooth regions range from elongate and irregular, to almost 
regular polygons. With strong contrast enhancement they show 
mottling on roughly the same horizontal scale of 15-30 m as seen 
in the waist, although they are less prominent than those markings 
within the waist. 

4.4. Jets 

Jets made visible by entrained dust and ice particles are promi- 
nent features of the region near the nucleus. Although the jets are 
not surface features, the approximate base locations are plotted in 
Fig. 3. These were mapped by interactive stereo viewing and plot- 
ting results on the projected shape model in single views of the sur- 
face. Because the contrasts of the jets change with viewing, and the 
background in most areas under the jets is variegated and rough, 
these positions, and even the number of jets, is subject to some 
interpretation. The basic occurrence pattern is reliable, and shows 
the association of jets with the rougher areas and mounds, and their 
avoidance of the smooth areas, as noted by A’Hearn et al. (201 1 ) and 
Syal et al. (2012). Accurate measurement of their dimensions is be- 
yond this work, but most are 3-5 pixels wide in the 7-m/pixel data. 
Detailed modeling of possible jet formation and relation to general 
surface evolution is reported in Syal et al. (2012). 

4.5. Other patterns and forms 

On the larger lobe, the region from ~260-320°E to ~20-70°N 
has somewhat rougher mound topography, more small spire-like 
mounds, and displays a somewhat higher albedo and bluer colors 
(Li et al., 2012) than the middle part of this lobe. The color differ- 
ences appear to arise from a greater exposure of water ice (Sun- 
shine et al., 2011; Li et al., 2012) and might relate to this area 
being near the morning terminator, the morphologic differences 
would not arise diurnally. In fact, there is a smaller region of some- 
what similar increased roughness in the region 30-80°E. Thus, 
although topographic type varies crudely along the long axis of this 
object, considerable variation exists longitudinally. 

5. Inferring a mean density 

5.1. Inferences from averaged quantities 

The distinctive, smooth profile of the waist was used by A’Hearn 
et al. (2011) and Richardson and Bowling (2012) to infer a mean 
density range for Hartley 2 using some assumptions. First, the 
smooth form was assumed to be at or close to an equipotential, 
presumably the result of depositional filling and smoothing a low 
area. Second, a simple 18-h rotation period was assumed about 
the maximum moment direction at the time of encounter. Third, 
the interior was assumed to be homogeneous in density. Minimiz- 
ing the variation of fractional potential energy over the area of the 



Fig. 7. Stereo views of parts of surface of Hartley 2. (a) Detail of the waist, contrast 
stretched to emphasize the variegated surface. North is to upper right, center is at 
approximately 55°N, 30°E. Panel width 430 m. Images mv3421 49303 J5000002 and 
mv342149307_5004051. (b) Region with darker, smoother areas and areas of small, 
brighter mounds on the larger, southern lobe. North is to upper right, centered at 
approximately 53°S, 31°E. Panel width 810 m. Images mv3421 49303 J5000002 and 
mv342 1 493 1 2_5004052. 


well-imaged waist suggested a mean density of 220 kg m -3 
(A’Hearn et al., 2011) and 140-520 kgm -3 (Richardson and Bowl- 
ing, 2012). Here we use the updated shape model, slightly different 
definitions of the waist area, and considerations of whether com- 
plications in the spin have practical effects to interpret the mean 
density. 

The most important additional concern is the comet’s complex 
and changing rotation (A’Hearn et al., 2011; Belton et al., 2012). 
First, the primary period was rapidly increasing at the time of 
encounter, approximately 0.1% per period; where the period was 
~16h only 60 days before encounter (Meech et al., 2011; Belton 
et al., 2012). Surface features are unlikely to evolve at meters/day 
(see Section 6.1), and thus the shape we measure and attempt to 
model would have formed over a considerable interval of time be- 
fore the flyby. Thus, the applicable primary period would be 16 h or 
possibly even less. Additionally, the complex spin imposes a roll 
about the long axis, which for our purposes may be thought of as 
having two effects: (1) The primary spin pole moves; in our map- 
ping the spin pole would move through different longitudes, vary- 
ing the rotational accelerations added to the gravitational 
accelerations. (2) The roll slightly reduces the relative gravitational 
topography along the longitude lines (along the cross profiles of 
the waist): For a 27.8 h roll period relative topography from the 
middle of the waist (latitude 40°N) to the northern edge (60°N) de- 




P.C. Thomas et al./Icarus 222 (2013) 550-558 


555 




Mean density kg nr 3 

Fig. 8. Variation of slope and potential as a function of mean density. Left: Average slope, degrees, as a function of mean density for the whole object, the entire waist (see 
Fig. 2), and the darker portion of the waist (Fig. 2). Right: variation of potential energy for same areas as in (a). 


creases by 1 m relative to a no-roll situation. Moving the spin pole 
position has some effect on the relative heights of points within the 
waist relative to the south pole because the waist is not perfectly 
circular along latitude lines (typically 25-30 m difference in radius 
from 60°E to 300°E along latitude lines in the waist). The north- 
south difference in dynamic heights along a particular longitude 
line varies by less than 1 m with the change in spin pole location. 
Thus, we can make comparisons based on a nominal spin pole loca- 
tion as the effects of the roll are effectively much smaller than is 
our ability to discriminate topography. We make these compari- 
sons of heights for different densities and spin rates. 

For each mean density we calculate the potential energy at the 
surface (Thomas, 1993), the effective acceleration, its angle to the 
surface normal (slope), and the dynamic height (Vanicek and Krak- 
iwsky, 1986) which uses an area-averaged acceleration to scale the 
potential energy values. We tested several areas for the variations 
in slope and potential: the whole object, the darker part of the 
waist (Fig. 3), the whole mapped waist, and the entire 360° longi- 
tude range between 10°N and 60°N. For these different areas of the 
surface we find the area-averaged slopes, and the area-averaged 
variance of the potential, the latter scaled for the average potential 
because different densities impose different absolute potentials. 

We made model predictions assuming primary rotation periods 
of 14, 16 and 18 h, and for bulk nuclear mean densities between 
120 and 800 kg m -3 . The results for the 16-h period, likely the 
most relevant, are shown in Fig. 8. The results for an 18-h period 
give minima at slightly lower mean densities; the results for a 
14-h period yield slightly higher mean densities. Results from cal- 
culation of average slope and variation of the potential energy 
(scaled for average values within an area) give similar, but not 
identical relations. Modeling the entire surface yields a distinct 
minimum for both average slope and potential variance (Table 2). 
These minima range from 265 kg m -3 to 300 kg m -3 . The waist re- 
gion gives a minimum average slope at a density of 340 kg m -3 , but 
yields no minimum using potential variance. The dark portion of 
the waist provides no minimum for any period or measure. The dif- 
ference in these results from those of Richardson and Bowling is 


Table 2 

Densities for minima of variance of slopes and of potential. 


Area 

Density from 
slope (kg m -3 ) 

Density from potential 
variance (kg m -3 ) 

All 

300 

265 

Waist 

340 

NM 

Dark 

NM 

NM 

Period = 16 h 



NM = no minimum 




chiefly in the area selected; results of whole-body measures in 
the two studies are essentially identical. Because these minima 
are judged significant based on process assumptions, we have no 
formal way of applying uncertainty ranges to these values in 
Fig. 8. From these plots alone we would infer a mean density be- 
tween 265 and 340 kg m -3 . 


5.2. Insights from topographic profiles 

We attempt to refine the density inferences by reference to 
topographic profiles shown in Fig. 9. These show north-south pro- 
files of dynamic heights at two longitudes for the 16-h rotation 
period. These profiles emphasize the variation of topography with- 
in the waist and the range of mean densities that come close to 
making part of the waist nearly flat, but which still leave noticeable 
slopes at other locales. 

The simplest approach is to see what mean density gives the 
waist area (generally on these plots between x-axis values of 
0.05 and 0.80 km) the flattest profiles; this value is ~230- 
240 kg m -3 . Values of 200-300 kg m -3 also appear consistent with 
a low, relatively flat waist region. These profiles also illustrate that 
global relief is minimal for densities near 250 kg m -3 . 

Do the possible geological interpretations of the waist warrant 
matching an equipotential? This smooth surface might represent 
an equipotential surface, that is, it might be an area of ponding 
of loose material that cannot support slopes and is all close to 
the same potential. It might represent a relatively low area that 
is accumulating materials but has not entirely been flattened. It 
is bounded on the south by a transition to rougher terrain, some- 
times gradually, at latitudes ranging from ~0 to ~30°N. The north- 
ern boundary is gradational over ~100 m at latitudes of ~56-60°N. 
The waist does not have a sharp boundary such as do “ponds” on 
Eros (Robinson et al., 2001) or the Muses Sea on Itokawa (Fujiwara 
et al., 2006). Although ponds are visibly distinct, there is a sugges- 
tion that some of the ponds on Eros are not completely “flat” (Rob- 
erts et al., 2012). The waist on Hartley 2 does not show flow fronts 
or other likely indicators of discrete flow units. Thus, while the 
images suggest low-slope processes, there is no clear geological 
restriction on how flat this region is. 

The other darker smooth areas in the large lobe might be ex- 
pected to either be flat, locally low, or perhaps show slopes indic- 
ative of a flow, especially for the elongate one at ~50°S, and 338- 
38°E (Fig. 3). These features occur in slightly lower regions than 
their surroundings, though the differences are largely <15 m. We 
note that the rotational light curves of H 2 0 and C0 2 (Besse et al., 
201 1 ) imply that the smaller lobe is much more active than the lar- 
ger lobe so that transfer of material from the smaller lobe to these 
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Fig. 9. Shape and dynamic heights along two longitudinal traverses. (A) 310°E, 
radius to surface. (B) Calculated topography along trace as a function of assumed 
mean density in kg m -3 , noted for three of the curves. Spin period of 16 h. (C) Same 
as (A) for 10°E. (D) Same as (B) for 10°E. The potential energies, or heights, are 
calculated relative to the south pole. Other points at the same rotational distance 
can all have the same rotational energy, which added to the same relative body 
mass potential results in curves recrossing at a point. 

regions is plausible. Slopes within most dark areas exceed 5° no 
matter what the assumed density or rotation period. 

5.3. Density summary 

The model results in Fig. 8 and Table 2 suggest a mean density 
of 265-340 kg m -3 . The two sample profiles in Fig. 9 suggest a 
mean density between 200 and 300 kg m -3 . The limitations of 
our assumptions mean that systematic effects might allow a 
broader range of mean densities, probably on the higher side. We 
adopt a range of 200-400 kg m -3 (300 ±100 kg m -3 ) as approxi- 
mate range using minimal significant figures. The upper value is 
not really constrained by Fig. 8; interpretation of Fig. 9 is the only 
upper bound possible with our approach, and we adopt 400 kg m -3 
because of Fig. 8 results and the possible applicability of a still 
shorter rotation period. These results are broadly consistent with 
those of Richardson and Bowling (2012) who obtain 140- 
520 kg m -3 . The inferred density of Hartley 2 is within the range 
of estimates for Tempel 1 of 400 + 600/-200 kg m“ 3 (Richardson 
et al., 2007). 

The density modeling uses only a homogeneous interior. If this 
is an assembly of two originally separate pieces, each might have a 
different mean density. Additionally, the deeper interior might 


have a different density from that of the surface regions. Investiga- 
tion of these parameter spaces is left for later study. The simplicity 
of the current model emphasizes that the inferred density of 200- 
400 kg m -3 probably is an optimistic rendering of the uncertainty. 
By far our greatest concern with this modeling apart from assump- 
tions of geologic processes, is that it is a static model of dynamic 
processes: the rotation and the surface accelerations have been 
changing, and the surface has been changing also. The observed 
shape is almost certainly an integrated effect of processes that 
has not reached, nor will it likely reach, an equilibrium form. 

5.4. Meaning of a very low density 

A density of 300 kg m -3 would require a high mean nucleus 
porosity, >70%, that could originate by inefficient packing of elon- 
gate or irregularly shaped particles and/or macroscopic voids. Par- 
ticle densities would be roughly 930 kg m -3 for H 2 0, 1500 kg m -3 
for C0 2 , ~800 kg m -3 for CO and 2500-3500 kg m -3 for refractory, 
rock-forming materials (e.g., silicates, metal sulfides, amorphous 
carbon). An upper limit of mean density of 520 kg m -3 would re- 
quire ~45% porosity for water ice particles, and >50% for modestly 
more dense particles such as C0 2 ice or silicates. The ejection of 
loose clumps of material that rapidly disaggregate (A’Hearn et al., 
2011) is consistent with a low-density assemblage. The different 
terrains on the surface indicate the operative processes sort mate- 
rials and/or have different stages of development. The different rel- 
ative amounts of C0 2 and H 2 0 (Besse et al., 2011; Feaga et al., 
2011) also indicate either sorting of materials, global inhomogene- 
ities, or different stages of processes in different regions of the co- 
met nucleus. 

6. Changing Hartley 

6.1. Mass loss by Hartley 2 

In order to estimate the total mass loss by Hartley 2 at a given 
apparition, we consider primarily the observations of the water re- 
lease as determined from Ly-oe measurements with the SWAN 
instrument on the Solar and Heliospheric Observatory (SOHO) 
(Combi et al., 201 1 a, b). These data cover the 1997 and 2010 appa- 
ritions and show clearly that the activity decreased by roughly a 
factor 3 from 1997 to 2010. The SNR is much better in 1997 so 
we use that apparition to integrate around the orbit. 

In the 1997 apparition, the water production peaked at 
log(Q.) = 28.55 molecules s -1 a few days after perihelion. We have 
carried out a numerical integration over the observed portion of 
the orbit using linear interpolation between each of the data points 
and we have also analytically integrated the power law obtained 
by assuming straight lines, eyeball fitted in log(Q) vs. log(At) (stee- 
per before perihelion than after; At is time since perihelion). The 
two approaches agree within 20%. We have extended the range 
in heliocentric distance slightly using two points based on OH 
measurements supplied by Knight and Schleicher (2012) at 
50 days pre-perihelion, log(Q.) = 27.64 s _1 , and 66 days post-peri- 
helion, log(Q) = 27.69 s _1 , both nearly an order of magnitude below 
the peak. To estimate the contribution further from perihelion, 
Meech et al. (2011) show from water-sublimation models, that 
the heliocentric light curve data (scattered light measurements 
from the dust) is consistent with the activity beginning ~450- 
500 days pre-perihelion near 4.3-4.2 AU. This is consistent with 
what was seen in previous apparitions by observers as compiled 
by Ferrin (2010). Post-perihelion activity continued through aph- 
elion, turning off between 200 and 400 days post-perihelion (near 
r = 5.6 AU). Although Meech et al. (2011) suggest that this contin- 
ued presence of grains in the coma is from C0 2 outgassing, we as- 
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sume here that the continuation through perihelion is due primar- 
ily to the rapidly increasing lifetime at large heliocentric distances 
of the icy grains that were released earlier and seen so prominently 
at this apparition (A’Hearn et al., 2011) and that the outgassing 
stops at 400 days. This yields a conservative lower limit on the to- 
tal mass loss. If the average production rate during the decline to 
zero is as high as 10 24 s _1 , roughly 1/3 the value suggested by 
the apparent magnitude at 250 days, it adds only about 2% at each 
end to the total mass loss. This computation yields a conservative 
lower limit of the total water release of 1.5 x 10 35 molecules or 
4.5 x 10 12 g. From measurements during the flyby we know that 
C0 2 is about 30% by mass of H 2 0. If we assume dust/gas = 1, the to- 
tal mass loss during the 1997 apparition was 1.2 x 10 13 g. During 
the 2010 apparition, therefore, it was roughly 4 x 10 12 g. These 
numbers are very sensitive to the dust/gas ratio, a poorly deter- 
mined parameter for all comets (due to the mass of dust being 
dominated by the largest grains in many comets), but half this va- 
lue (the ices) is a solid lower limit and the value could range up to 2 
or 3 times larger than the value above. 

Using the surface area from Table 1, the mass loss in 1997 cor- 
responds to ~2300 kg m“ 2 . For a bulk density of 300 kg m“ 3 , this 
value corresponds to a loss of 7.7 m of material everywhere on 
the surface. In terms of the upper range of the density from Rich- 
ardson and Bowling (2012) (and a more typical density for other, 
better determined nuclei) a mean density of 520 kg m -3 would 
suggest loss of 4.4 m per orbit. Using the volume from Table 1, 
and a density of 300 kg m -3 , the mass loss during the 1997 appa- 
rition was roughly 5% of the total mass of the nucleus and thus 
about 2% of the total mass during the 2010 apparition. An upper 
density value of 520 kg m -3 would suggest a loss nearer 3% of 
the nuclear mass in 1997 and 1% in 2010. Even allowing for a factor 
of a few uncertainty due to the unknown dust/gas mass ratio, this 
mass loss, both as an average depth and as a fraction of the mass of 
the nucleus, is far higher than that found for the much larger nu- 
cleus of 9P/Tempel 1 (~l/3 m/orbit or ~3 x 10“ 4 mass/orbit; Tho- 
mas et al., 2012). 

The above calculations are average amounts over the entire sur- 
face. However, all the evidence from the flyby (A’Hearn et al., 201 1 ) 
suggests that much of the material is lost from sub-surface out- 
flows that are localized, so collapse features might be the expected 
expression one would see in images from different apparitions, 
perhaps analogous to those seen near the smooth flow on 9P/Tem- 
pel 1 (Schultz et al., 2011). However, collapse forms are not obvi- 
ous in the images. Perhaps some of the dark areas might be 
termed sags (Section 5.2) related to loss of material from extended 
sub-surface regions. 

6.2. How is the surface shaped? 

The somewhat more elongated average shape of cometary nu- 
clei compared to small asteroids has been attributed to effects of 
surface mass loss by sublimation (Jewitt et al., 2003). In detail, 
however, sublimation lowering of surfaces generates positive feed- 
back effects (such as formation on Earth of suncups, penitents, or 
inverted topography; Lliboutry (1954), and discussed in the case 
of comet nuclei by Malin and Zimbelman (1986) that leave a rough 
and variegated surface quite unlike what is visible on Hartley 2. 
Our close-up views of comets are limited and in situ observations 
nonexistent, so we must consider general principles and make the 
best available comparisons to infer how this object evolves. 

The bi-lobed shape suggests possible formation from two 
pieces. The near axial symmetry suggests effects related to spin 
and gravity rather than to random events such as impact cratering. 
As noted above this comet does not display the erosional and depo- 
sitional forms found on Tempel 1 (and likely on Wild 2 and Borrel- 
ly; Soderblom et al., 2002; Brownlee et al., 2004) that contribute to 


a complicated shape that is not symmetric and which has consid- 
erable topography relative to gravity. Indeed, even though dis- 
persal of material in jets, carrying particles that continue to 
fragment (A’Hearn et al., 2011) is obvious, the vent geometry, if 
there are discrete vents, is elusive, and any flow forms likely vastly 
different from those on Tempel 1 (Belton and Melosh, 2009). There 
are some areas, the darker patches (Section 4.3) that suggest depo- 
sition, but which are not gravitationally confined materials such as 
the “ponds” on Eros, because they have some gravitational slopes. 
Even the waist area has roughness on the scale of meters vertically 
and regional slopes of a few degrees for any mean density, and thus 
is also not a simple deposit, even though its shape strongly sug- 
gests some mechanism of regional smoothing. Thus we must look 
for more indirect routes of molding the surface of this comet. 

As noted in Section 6.1 the mass flux from this comet could 
amount to several meters of surface lowering per orbit, and given 
the distributed jets and lack of large vent morphology, the possibil- 
ity arises of subsurface cavities supplying material and then being 
subject to collapse. The size and depth of such cavities would pre- 
sumably not greatly exceed an orbital skin depth, thus would need 
to be numerous, nearly ubiquitous, and only marginally resolvable. 
At 7 m/pixel the surface shows no indication of dominance of col- 
lapse morphology. A trough <15 m deep exists at high southern lat- 
itudes (Fig. 1 bottom lit area; Fig. 7 near bottom), but much of this 
topography may be constructional clustering of the bright knobs. 
More intriguing are the smooth, slightly sloping, darker areas that 
appear to be either parts of bowls or of tilted segments, albeit near 
the limits of our measurement precision (Fig. 7). These depressions 
might be remnants of collapse of connected large voids partially 
filled and smoothed by material formerly above the voids (“crust”) 
or partially covered by some materials mobilized near the surface 
by other vents. 

The presence of mounds may suggest eruptive and constructive 
mechanisms at work, and the ejection of particles that subse- 
quently supply much of the mass loss by sublimation (A’Hearn 
et al., 2011; Harmon et al., 2011; Kelley et al., 2012) indicates an 
environment rich with material that could form brighter local 
deposits. There is the possibility that the average surface is losing 
material while some of the mounds undergo deposition and up- 
ward growth. The height of the mounds (generally <30 m; Sec- 
tion 4.1) might be the net result of loss and growth before the 
mounds collapse. This height might reflect the strength of the 
mounds, which can be crudely estimated based on pgh, where 
p = 300 kg m -3 , g = 4x1 0 -5 ms -2 , and h = 3 x 10 m, assuming a 
density at the low end of our mean values. This approximation 
gives a value of 0.4 Pa. Such a weak material would be consistent 
with local accumulation of what might be close to snowflakes. Gi- 
ven the small relief elsewhere on the comet nucleus, strengths esti- 
mated using these assumptions for materials anywhere on this 
object will be very low. 

The high rate of mass ejection does suggest there should be 
some small, low-velocity fraction that may not initially escape. 
This material might be distributed widely over the nucleus, and 
could possibly flow with the tenuous gasses evolving from ice par- 
ticles into the gravitationally lower regions. Such material would 
have essentially zero strength and be consistent with minimal 
topography on a very low-gravity object. 


7. Summary 

Hartley 2 is distinguished by a nearly axially symmetric, bi- 
lobed shape. An extremely active comet partly by virtue of expel- 
ling solid water-ice particles, its surface may suffer average loss 
of a few meters per orbit. Its surface features, chiefly mounds 
<30 m in dimension, relatively darker smooth areas, and a waist 
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with local topography <2 m do not easily compare to obvious ero- 
sional and depositional forms such as the pits and smooth deposits 
on Tempel 1. The smooth shape, and especially the waist, suggest 
processes whose net effects reduce topography on horizontal 
scales of >100 m. The shape combined with assumptions of 
smoothing processes suggests a mean density of 200-400 kg m -3 . 
This mean density implies loss of >1% of the nuclear mass per orbit. 
The axial near symmetry suggests no significant role for random 
events such as craters but does suggest interacting roles for insola- 
tion, spin, and gravity. 
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